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Abstract: Fluid dynamics corresponds to the dynamics of a substance in the long 
wavelength limit. Writing down all terms in a gradient (long wavelength) expansion 
up to second order for a relativistic system at vanishing charge density, one obtains 
the most general (causal) equations of motion for a fluid in the presence of shear 
and bulk viscosity, as well as the structure of the non-equilibrium entropy current. 
Requiring positivity of the divergence of the non-equilibrium entropy current relates 
some of its coefficients to those entering the equations of motion. I comment on possible 
applications of these results for conformal and non-conformal fluids. 
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1. Introduction 

The theory of fluid dynamics has a long history, starting with its inception by L. Euler 
in 1755 and generalized to viscous fluids by the works of C. Navier |2| and G. Stokes 
in the 19th century. Besides being necessary to realistically model the behavior of 
real fluids in many situations, viscosity is known to be essential for the presence of 
laminar (or smooth) flows since it dampens the turbulent instability generally inherent 
to ideal (non- viscous) fluids f|§26. As can be understood from the work by Navier and 
Stokes, viscous effects correspond to (first order) gradients of the equilibrium properties 
of the system (temperature, fluid velocity, etc.). 
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For relativistic systems, additional complications arise: a first order gradient ex- 
pansion a la Navier and Stokes leads to a set of fluid dynamics equations that allow 
faster-than-light signal propagation, violating causality ||. In the late 20th century, 
motivated by questions in general relativity, Miiller [||, Israel and Stewart [0] showed 
that by including second order gradient terms, the resulting fluid dynamics equations 
could be made causal. In the last decade - driven by nuclear physics experiments 
on relativistic heavy-ion collisions || [| [HJ 11 — a program to determine the precise 



form of the equations for a relativistic viscous fluid was initiated [[T|, [13|, [14], |15| . This 
program has not been finished, and the present work is meant to provide another step 
towards its completion. 

It should be noted that the question of causality is intimately linked to the property 
of hyperbolicity of the fluid dynamic theory, i.e. the question of whether a well-defined 
initial value formulation can be given. While for small perturbations around equilib- 
rium hyperbolicity and causality can be shown for various second order theories, it 
seems that in the case of strong non-linear out of equilibrium situations hyperbolicity 
has been rigorously investigated only in so-called divergence type theories [TBI, e.g. 



Miiller's theory [[J. In particular, the causality properties of the theory by Israel and 
Stewart far out of equilibrium remain unknown. This work is based on the more recent 
interpretation of fluid dynamics as an effective theory of the long-wavelength modes of 
the system, by construction precluding any application to systems that are far from 
equilibrium. While the effective theory specified below will share with other second 
theories the property of causality (and hyperbolicity) for small perturbations around 
equilibrium, a proof of even hyperbolicity for strong perturbations is beyond the scope 
of this work. 

This article is organized as follows: in section |2], all relevant gradient structures 
up to second order are listed for a fluid at zero charge density, and the concept of 
a conformal fluid is introduced. In section || the most general form of the energy- 
momentum tensor for non-conformal fluids is presented, which fixes the equations of 
motion for a relativistic viscous fluid. Section f| deals with the most general form 
of the entropy current for a relativistic fluid out of equilibrium (again, at vanishing 
charge density), and its divergence. In section |5|, I tried to collect all current knowledge 
about the coefficients multiplying second order gradient structures. Section |] contains 
a discussion of the results and the conclusions. 

2. Setup 

Let us consider matter described by a relativistic quantum-field theory at vanishing 
charge density (zero chemical potential). In equilibrium, the long- wavelength dynamics 
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of this system can be described by one scalar, one vector, and one tensor: 

which are the fundamental fluid dynamic variables energy density e, fluid four velocity 
and metric tensor g^ u . The quantum field theory is supposed to furnish a relation 
between the pressure P and e via the equation of state, P = P(e), giving rise to the 
speed of sound c s = a/ dPj de in the fluid. 

Using these fundamental degrees of freedom, one can write down the energy- 
momentum tensor of the fluid (see, e.g. Q §133) , 

T% = e uV + PA"" , A"" = + uV , (2.1) 

where the subscript "eq" denotes equilibrium quantities containing no gradients (zeroth 
order). Here and in the following the metric sign convention (— , +, +, +) is used, and 
all the calculations will be done in four space-time dimensions (though a generalization 
to other number of dimensions should be straightforward). Without the presence of a 
source, the energy-momentum tensor is covariantly conserved, V^T"" = 0, where V M 
denotes the geometric covariant derivative. Projection of this equation leads to the 
well-known equations of motion for an ideal relativistic fluid, 

UvV^Tg = -De - (e + P)V • u = , 

A^V M T^ = (e + P)Du a + VIP = , (2.2) 
where the new notations 

D = u»V^ , = A"" V„ 

for the comoving time and space derivative were introduced. Using the basic thermo- 
dynamic relations Se = TSs, e + P = sT for the equilibrium entropy density s and 
temperature T, the equations of motion for the ideal fluid become 

Dlns = -V-u, Du a = -c\ In s . (2.3) 

This implies that not all first order gradients of the fundamental degrees of freedom are 
independent: time derivatives may (up to higher order gradient corrections) always be 
recast as space derivatives. Therefore, to first order in gradients, the only independent 
structures one can write down are In s and V^w^ (no coordinate-invariant first 
order gradient of the metric tensor exists). For later convenience, gradients are sorted 
into three classes: scalars, vectors orthogonal to w M , and symmetric traceless tensors 
orthogonal to v?. To first order in gradients, there is one of each class: 

(V • u) , V^lns, a aP = + V^u a - - A Q/3 (V ■ u) . 
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To second order in gradients, the independent structures are V^Va Ins, V^V^u^, 
V^lnsV^lns, V^m^V^m^, V^w^V^lns and the Riemann tensor (cf. [p~7|] §3.4) 

pA __ o -pA o -p\ I -pK pA -pA 

flCTV — U V L 11V U V L ficr I" 1 /iZ/ 1 KCT 1 [1(7 KU > 

where P^, = (<9 M 5f piy + - (9^) . Using the Ricci tensor = and 
Ricci scalar R = R^, one can build 7 independent scalars to second order in gradients, 

ViVJlns, V^lnsV^lns, cx Q/ ,cx^ , (V ■ m) 2 , fl, 

(2.4) 

where Q^ u = ~ (V^-u„ — V^rw^) is the fluid vorticity which was used in the decomposi- 
tion 

v ij u u = -a^ u + i H — A^ u . 

Furthermore, one can build 6 independent vectors orthogonal to the fluid velocity: 

V> a ^, V^T", a Q ^lns, (V-u)V^lns, A a ^R a p , (2.5) 

where one should note that A^u^R^ contains - and therefore is used instead of - the 
vector (V • u) (see appendix |B|). 

Finally, there are 8 independent symmetric traceless tensors orthogonal to the fluid 
velocity. These have been found already in Ref. jT3| for the case of conformal fluids 
(see below for the definition of a conformal fluid). For a general relativistic fluid at 
vanishing charge density, the 8 independent tensors are given by 

V^lnsVflns, V^Vflns, (V • u) , a K V >A , 

a <fl x Sl u>x , Q <fl x Q u>x , u a u p R a< » v> P , R<^> , (2.6) 

where for a second rank tensor A^ u 

(A» v ) = A <iiV> = -A>* a A ul3 {A af} + Ap a ) - -A^A afS . 

2 3 

2.1 Conformal Fluids 

For certain situations it is advantageous to consider the simplified case of a fluid without 
bulk viscosity. Since the bulk viscosity coefficient £ is related to the conformal anomaly 
via a Kubo relation, a prime example for such a case is a system that exhibits 
conformal symmetry, or covariance under local Weyl rescalings of the metric: 

g^^e-^g^, (2.7) 

where w = w(x^) is the local scale factor. I will refer to the fluid description of a 
system that obeys conformal covariance simply as "conformal fluid" . An example for a 
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conformal quantum field theory is the M = 4 supersymmetric Yang-Mills theory (SYM). 
Note that in curved space the Weyl anomaly |18| would in general break conformal 
symmetry, since, e.g., oc R 2 . However, this breaking occurs only at fourth order in 
gradients, suggesting that for a fluid dynamic expansion up to third order in gradients 
the Weyl anomaly can be ignored [[14] . 

Conformal covariance implies that objects transform homogeneously under Weyl 
rescalings, meaning that gradients of the scale factor w have to cancel. To lowest order, 



consistency requires [[14 



s e 3w s , e w u" . 

To first order in gradients, there are no scalars or vectors 1 but one symmetric tensor 
that transforms homogeneously under Weyl rescalings: a^ u (see appendix 0). To second 
order in gradients, there are three conformal scalars, 

Si = , s 2 = , 

= c 2 s In 8 + f In sVt In s - \u a u p R a ? - \R + \ (V • uf , (2.8) 
two conformal vectors orthogonal to it M , 

VI = Via a » + 2c 2 s a a ^Vi In s - , V£ = + u^ a ^ , (2.9) 

and five conformal symmetric traceless tensors orthogonal to , 

Qi w = R< ^> _ c 2 (2V< M Vf In s + (V • u) - 2c 2 V< M In sVf In s) , 

Of = R<^> - 2u a upR a< ^» 3 , 

Of = a <\a v>x , Of = a K \Q u>x , Of = fi^fi"^ . (2.10) 

These will be the building blocks of the energy-momentum tensor and entropy current 
for conformal fluids. 



2.2 Non-Conformal Fluids 

For more general fluids that do not obey conformal invariance, all possible gradients 
can contribute. In particular, to first order in gradients there are one scalar and vector, 

(V-tt), V M lns, (2.11) 

in addition to the tensor already found for the conformal case. 

1 For non-vanishing charge density, parity-breaking effects may occur in the fluid. In this case, aiso 
terms such as e^"^ u v Q, a p are allowed. 
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At second order, there are four additional scalars, 



S A = (V-uf, S 5 = R, <S 6 = V^lnsV^lns, S 7 = u a u f, R a p, (2.12) 

four additional vectors, 

V 3 M = A^u a R a p , V 4 M = a^Vi In s , V 5 M = Q^Vi In s , V 6 M = ( V • u) In s , 

(2.13) 

and three additional tensors, 

0£" = u a u p R a<flu> P , Ot; u = ^-^a^, 0%" = V^lnsVf Ins. (2.14) 

3. Viscous Fluid Dynamics: Energy-Momentum Tensor 

The energy-momentum tensor for a relativistic viscous fluid can be written as 

± non— eq - 1 eq ~r 11 i 

where the viscous stress tensor IP 1 * contains correction terms to the ideal energy- 
momentum tensor due to shear and bulk viscosity. At this point, I recall that for 
vanishing charge density 2 the only useful way to define the fluid velocity is the 
Landau-Lifshitz condition 

u,T^ = eu\ 

which implies u fl U flL ' = 0. The viscous stress tensor is customarily separated into a 
traceless part (tt^) and a part with non- vanishing trace (II), 

W u = tt^ + A^IT . 

For conformal fluids, where = 0, the trace part vanishes identically and the 
structure of the traceless part tt^ in the fluid dynamic (small gradient) regime is 
generated by a^ u and the tensors found in Eq. fl2.10| ) (cf. Ref. JI4"|). For general fluids, 
also the tensors in Eq. Q2.14| ) contribute. Using instead of 0±* the expression 3 



9 — V 2 1 rl r 2 

I Da^) + -^a^ (V • u) ~ OT - Of - -Of + 2CT - In s Vf In s , 

v / 3 1 2 2 A 5 dlns ± ^ 

(3.1) 



2 For non-vanishing charge density p, the presence of the charge current j M offers other choices to 
define u M , such as u^j^ = p. In the Landau-Lifshitz frame, since is defined via the rest-frame of 
the energy-density, heat diffusion does not exist; but, since in this frame generally u^j^ =/= p. there 
can be charge diffusion. In the so-called Eckart frame u^j^ — p, and there is heat diffusion instead of 
charge diffusion, indicating that these concepts are different manifestations of the same phenomenon. 

3 Note the wrong sign of the O5 term in Ref. juj that has been corrected here. 
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which is accurate to second order in gradients and amounts to a particular resummation 
of higher order terms, one finds 



+A 1 a < V" >A + X 2 a < \n v>x + X 3 n <tA x Q u>x 

V • 11 

+K*2u a iipR a<flu> P + 7]rl + A 4 V <M In sV > In s . (3.2) 

3 

The coefficient of the first order gradient term is the familiar shear viscosity coef- 
ficient rj. The coefficients 7V, r*, k, k*, Ai, A2, A3, A4 are "second order" transport coef- 
ficients, three of which (t*,k*,A4) must be identically zero for conformal fluids since 
they multiply structures that do not transform homogeneously under Weyl rescalings. 

The expression for the trace part IT in the fluid dynamic regime contains all scalars 
up to second order in gradients, which are given by Eq. fl2.8| ) and Eq. ( |2.12| ). Using 
instead of S 3 the expression 

D (V.«) * -\s 1+ s 2 -s 3 - l -s, - is 5 + - s.-**,, 

one finds 

n = -C (V • u) + (r u D (V • u) + £i<t^<v + 6 (V • uf 

+£ 3 C^CV + U^t In sVl In s + + i & u a u?R aP . (3.3) 

Here ( is the familiar bulk viscosity coefficient and rn, £1, £2, £3, £4, £5, £e are second 
order transport coefficients for non-conformal fluids. As will be shown below, at least 
two of these second order coefficients (e.g., £5,£6) are completely specified in terms of 
other transport coefficients. The equations (|3 . 2|) , (|5T3|) give the most general structure 
for the energy-momentum tensor of a relativistic viscous fluid at zero charge density 
up to second order in gradients. 

Similar to the case for ideal fluids, the equations of motion for a viscous fluid are 
given by V^T^ = 0. In particular, using again the basic thermodynamic relations one 
finds for uJW ^T^ u = the result 



Ds-i ,.(V.«) = ~7r^^-^n(V- ( /) . (3.1) 



which will be useful in the following. 
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3.1 Dispersion Relations and Kubo Formula 

Considering small perturbations 5s, 5u^ around an equilibrium configuration in flat 
space with s = const, and = (1,0), one obtains dispersion relations for the collective 
modes along (L) and perpendicular (T) to the perturbation. The analysis is straight- 
forward (see for example Ref. |19| §11. C) and in the fluid dynamic regime (u,k <C 1) 



leads to 

k 3 

oj L (k) = ±kc s - ik 2 T =f — 

2c, 



r 2 - 2c 2 1 - 



2 T]T n 1 (t h 



3s + P 2s + P 



+ 0{k 4 



u T (k) = + 0{k') , (3.5) 

where it is recalled that transverse and longitudinal perturbations correspond to the 
shear and sound mode, respectively, and the "sound attenuation length" is given by 

V 1 C 



3s + P 2s + P 

The equations of motion for a viscous fluid are found to be causal if the maximal 
propagation speeds at high wavenumber k for the sound and shear mode are smaller 
than the speed of light (cf. |19[ §11). Using the dispersion relations, one finds for the 
propagation speeds in natural units (J% = c = k B = 1): 



lim dUT{k) = I V lim dUL{k) = J<? + - 5 + < 

fc^oo dk ~ \jr n {s + P)' k^oo dk y s 3T n {s + P) r n {s + P) ' 

(3.6) 

The values are fixed by the first and second order transport coefficients, or the 
properties of system in the fluid dynamic (long wavelength) limit, while causality con- 
cerns the property of the system for short wavelength perturbations. There is a priori 
no reason why lim^oo UJ / J ' L should be less than unity 4 . 

The second order transport coefficients are related to thermal correlators via general 
Kubo formulas. As an example, let us calculate the retarded correlator G^ ,xy for the 
energy momentum tensor component T xy '. Considering a metric perturbation 5g iW with 



4 One could still use the second order equations of motion as a phenomenological model of a relativis- 
tic viscous fluid if the above conditions were violated for a particular quantum field theory: adjusting 
the second order transport coefficients 7V , rn "by hand" to repair causality, the resulting model would 
still treat the first order gradients correctly. In a way, something similar is done for numerical simula- 
tions of ideal fluids: to dampen the turbulent instability inherent to ideal fluid dynamics, one has to 
introduce "numerical viscosity" (first order gradient correction terms) to obtain stable evolution for 
the ideal (zeroth order gradient) fluid. Here second order gradients are needed to obtain stable and 
causal evolution of a viscous (first order gradient) fluid. 
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only the non- vanishing component 5g xy (t, z), the fluid stays at rest and in equilibrium: 
e = const, = (1,0), since this corresponds to a tensor perturbation. The correlator 
may then be found from the linear response of T xy to this perturbation (cf. HT4| ): 

G xy, X y = p _ irj u + ^ ^ _ « + ^ ^2 _ « fc 2 + fe 3) _ 

This implies that once tv is known, both k and n* can be found by calculating this 
correlator for a quantum field theory in the appropriate limit. 



4. Non-Equilibrium Entropy 

For a system in equilibrium the entropy current is simply given by the product of 
entropy density and fluid velocity, S£ q = sw M . For systems out of equilibrium it may be 
that the entropy current gets modified 5 . It is known from kinetic theory that at zero 
chemical potential (in the absence of heat/charge diffusion), the first correction to the 
equilibrium entropy current must be of second order in gradients [PJ. In the past, the 
form of the non-equilibrium entropy current has often be postulated [0 . More recently, 
a more fundamental approach has been advocated |2l], |22] that calls for all structures 
in a gradient expansion to be allowed. I will follow this approach here, recovering and 



extending some of the results in Ref. [22 . 
4.1 Conformal Fluids 

For conformal fluids, the entropy current must be built out of elements that are invariant 
under conformal transformations, which are the three scalars (|2.8|) and two vectors ( [2.91) : 



+Bi Qvf + ^<S^ + B 2 (V£ - , (4.1) 

where the five coefficients ^1,2,3 and are (mass dimension one) functions of entropy 
only and the combinations and prefactors have been chosen such as to facilitate compar- 



ison to Ref. 1 12 2| I . For conformal fluids in three dimensions one has ct = | and II = 0, and 



the absence of a second dimensionful scale leads to the relation A4, ~ s 1 / 3 . Accord- 
ing to Boltzmann's H-theorem, entropy is never allowed to decrease, so the divergence 
of the non-equilibrium entropy current should obey 

xi ^ > n 

v M non— eq — u • 



3 I am not aware of any experimental verification of this hypothesis. 
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The divergence of the entropy current is a physical observable, and as such should 
transform homogeneously under Weyl rescalings. Explicitly, one can convince oneself 
that this is the case by writing 

1 P 4 ™ 

Taking the covariant derivative of Eq. ( f4.1| ), the result for the equilibrium part 
V /t (su^) can be read off from Eq. ( |3.4| ). A somewhat more lengthy calculation (see 
appendix [FJ for some useful identities) gives 

V^-eq = \K^W V (-2A 3 + B x ) + In a (-2A 3 + B x ) 



+0. 



fJtV 



^ + RT (— + A 3 ) + u a u p R*<^ ( ^p- +A 1 + B 1 - 2A 3 



-\o» x o» x r Xl ~^ +A 1 + B 1 - 2A 3 j + IvTvr In s - A x - 2/ 
+W a Qr (^-h±pil + Al _ 2A 2 - 2A 3 + B^j (4.2) 

where the conformal fluid property r* = k* = A4 = was used. Note that B 2 completely 
drops out the the divergence of the entropy current because the covariant derivative 
of the relevant term is vanishing. Positivity of V A1 S , Q 0n _ eq is usually guaranteed by the 
term a^a^, which is of second order in gradients and hence generally much larger 
than all the other (third order gradient) terms. However, there is the possibility that 
a^y itself is accidentally small, e.g., when considering a fluid velocity field that has very 
little shear motion in it. In this case it can happen that third order gradients dominate 
the entropy production, and hence their coefficients must be such that entropy never 
decreases. This immediately implies 

B 1 = 2A 3} (4.3) 

since the first two terms of Eq. (|4.2| ) could otherwise lead to negative entropy production 
(this relation was already pointed out in Ref. p2| ). Now by the same logic, one would 
also expect the coefficients of the term a^ u R^ u to vanish, leading to 

A3 = £ ■ (4.4) 
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However, there is a possible loophole in the argument leading to Eq . ( |4.4|) : it could be 
that the term a^ u R IJjU combines with other terms of second and fourth order in gradients 
to form a full square, e.g., 



a 



^-o^ + R" U (a 3 - — 
L2T V s 2T 



■—(a 3 -—) 2 r 11u r^ 

2r] V s 27V M 



2T 



a 



2t) R ^ 



a 



. (4.5) 



In this scenario (pointed out in Ref. ||22||), the coefficient A 3 could be arbitrary and 
positivity of the entropy current would still be guaranteed if the offending first term on 



the r.h.s of Eq. ( }4.5|) is offset by another term of fourth order in gradients. There are 
three possible sources of fourth order gradient terms in Eq. (f4.2|) : third order gradients 



in 7i tiu times stemming from V M (sit^), first order gradients in ^ v times third order 



gradients stemming from viscous corrections to Eq. (|2.3| ), and finally derivatives of third 
order gradients in the non-equilibrium entropy current S^ on _ eq . First note that the term 
R^yR^ cannot be offset by a fourth order gradient term from V M (sw M ) because — as 
Eq. 
Ru 



shows — the latter always involves a^ u which is not expressible in terms of 

73|) up to second order in gradients are 



Secondly, viscous corrections to Eq. 
found to be 



D In s ~ — ( V • u) + 7] 
Du a ~ -cl V7 In s + 



fW 



2sT ' 
sT 



/JO 



sT 



U a T) 



fJ,U 



2sT 



(4.6) 



which shows that — again — the term R PjV R^ v cannot be offset by these contributions. 
Finally, all third order gradient terms contributing to the entropy current may be found 
by the same principle as in section |[ Explicitly, one finds as scalars: 

• cr^v times all symmetric second order tensors (e.g., a^O^), 

• V^lns times all second order vectors (e.g., V^lns Vf), 

• (V • u) times all second order scalars (e.g., (V • u) Si), 

• acting on all second order vectors (e.g., V^Vf), 
and as vectors: 

• V^lns times all second order tensors (e.g., Ins O^), 

• (V • u) , a^ v , VL^y times all second order vectors (e.g., (V • u) Vf , a^ a V"), 
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• V^lns times all second order scalars (e.g., V^lnsSi), 

• acting on all second order scalars (e.g., V^<S>i), 

• acting on all second order tensors (e.g., V^C"^) . 

One can readily convince oneself that derivatives of third order vectors in the entropy 
current will not be able to generate a term such as R^ V R^ U . For scalars, the only 
potentially dangerous term is the time derivative of a^ v , because contains R^ v 
and Da^y could potentially contain another R^. But the explicit expression fl3.1|) for 
Da^ u shows that this is not the case, hence no R^ U R^ V term in V^S^ on _ eq is generated 
from any of the possible sources. 

As a consequence, the scenario Q4.5|) cannot be possibly realized for every A 3 . Only 
the condition ( |4.4|) ensures positivity of the non-equilibrium entropy. 

A similar (but slightly more complicated) argument can be made for the terms 
u a u / 3R a<flu> ^u' y u s R^ < ^ u> s and Q <lia Q u> a Q < ^Q u J! : for these terms to be canceled or 
made positive definite, one has to add the two third-order gradient terms x^o^ {0^ u — O^) 

iiv L <J ^non-eq 



and x 2 (Ti tl , 20^ u to S^ on with the requirements 



To cancel the cross-term u ol upR a<fl ' />l3 Q <tJ, ' y Q u > , however, implies the condition 

The resulting inequality relation can only be fulfilled if 

2k + A 3 . . 

Unfortunately, I did not find an argument that would fix the value of A\ in 
terms of second order transport coefficients. Nevertheless, it seems that a bound on 
the value of A\ could be found by extending 71^ to third order in gradients, since 
xxo^D {Oi V — Orf) involves D {u a u/3R a<>lu>li ^j which seemingly has to be canceled 
exactly by a corresponding expression in tt^ v , fixing x\. 

To conclude, requiring positivity of the divergence of the non-equilibrium entropy 
current fixes three of the five possible coefficients (A2, A3, B{). The divergence then 
takes the form 

V7 cm ^ „ „iiv 1 K ~ 2Al a uX 



'A 1 K-T]T n 



(Da^) + -a^ (V ■ u) 
3 



(4- 
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which fixes the rate of entropy production slightly out of equilibrium for any conformal 
fluid up to one unknown parameter (Ai). Note that for any theory with k 7^ 2Ai, the 
divergence of the non-equilibrium entropy current must differ from the simple expec- 
tation = ^<V<7^. 

4.2 Non-Conformal Fluids 

For non-conformal fluids, the non-equilibrium entropy current may be built out of all 
gradient structures. Keeping the form of the conformal entropy current intact, up to 
second order in gradients this leads to 

Snon-eq = su<* + + A 2 S 2 u» + A 3 (4S3 - ~<Si + 2S 2 ^j u" + A A S^ + A 5 S 5 U» 

+A 6 S 6 u» + B x Qvf + ^SA + B 2 (V£ - u»S 2 ) + B 3 (v» - ^VfJ + B 4 V£ 
+B 5 Vg + B*V$ + 4 1 ( ~ 2 S^ + \s 3 u» - ^- 2 S^ + - 2 S^ + 2(? s S 6 u» + 3Vf 

+6V g _ 6 y£ + (4 - Qc 2 M - ^c 2 X - \vi\ , (4.9) 

where the cumbersome expression multiplying A 7 leads to a particularly simple diver- 
gence (see appendix |FJ). 

Since for non-conformal fluids c 2 may be a function of the entropy density, its 
derivative c' s = ^jj^f is in general non-vanishing. Similarly, the coefficient functions 
Ai, Bi now may contain logarithms or any power of the entropy density, so that in 
general A' { = d Jf ^ ^ | (and the same for Bi). Using Eq. (|3.4|) , the divergence of the 
non-equilibrium entropy current can be calculated, but the resulting expression (|C.2|) 
is suitably lengthy to be unenlightening except maybe for the expert reader, so it has 
been relegated to the appendix. Nevertheless, the principle of positivity of the entropy 
current divergence singles out conditions for the Ai , Bi as was the case for conformal 
fluids. One immediately finds 

A 5 = 0, B l = 2A 3 , B 3 = 2(l-3c 2 )A 3 , B 4 = -A 3 (A 3 - 3c 2 (l - 2c 2 )) , (4.10) 
B 5 = B 2 (B' 2 - c 2 s ) + (B> 3 + c 2 )B 3 , B 6 = 2A 6 - A 7 (A 7 ~ <2) + K 1 ~ ?>c 2 s )A 3 B' 3 . 



In addition, one can again go through the arguments as in section |4.1| to show that 
o^W"" and (V • u) R cannot form full squares because no other contribution to the 
divergence of the entropy current could offset their negative contribution to fourth 
order in derivatives. This leads to the conditions 

K A 3 T , 

A 3 = — , 6 = -!_(34-l). 
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Finally, requiring that the cross-terms w Q! w i a.R Q<A "' >/3 f2 <At7 f2 J ,> and u a u^R a ^Q flu Q tlu can- 
cel leads to 



6/c(l — 2ci) — 2k* + A 3 + ^2/ 9 ,/n / 9n , / 9 /. 

^2 = * ^ , = -^(l-Qc 2 s +3A' 2 )+2(l-3cl)A 3 (2c 2 s -B' 3j 



Contrary to the situation for conformal fluids, fixing 8 of the 13 coefficient functions 
in the entropy current still does not seem to lead to a simple form for V^S^n.^. It 
may be possible that a more detailed analysis of the effect of third order gradients S 11 , 
7i fiu and II could result in further conditions on either the remaining five coefficient 
functions or some of the second order transport coefficients r*, k*, A4, £1, £2, £3, £4- 



5. Second Order Transport Coefficients: Known Results 

The most general structure for the viscous energy-momentum tensor up to second order 
in gradients has been established in section [| For conformal fluids at vanishing charge 
density, the energy-momentum tensor — and hence the equations of motion for viscous 
fluid dynamics — depend on 7 dimensionless numbers: the speed of sound and the 
transport coefficients. To first order in gradients, the only transport coefficient is the 
shear viscosity 77, and to second order there are in general five additional coefficients: 
TV, k, Ai, A 2 , A 3 . All of these transport coefficients seem to be independent, but so 
far A3 = has been found for all examples of quantum field theories where it has 
been calculated, perhaps suggesting that there is another symmetry that has not been 
exploited yet. 

For non-conformal fluids, the energy-momentum tensor can no longer be expressed 
in terms of dimensionless numbers because of the presence of an additional scale. Com- 
binations involving transport coefficients that have vanishing mass dimension then are 
functions of, e.g., the entropy density. Also, there are 9 new independent transport 
coefficients in addition to the 6 for conformal fluids: to first order in gradients the bulk 
viscosity ( and to second order rn, r*, k*, A4, £1, £2, £3, £4, while £5,^6 are specified in 
terms of the other transport coefficients and their derivatives. 

Many studies have been published on the value of the first order transport coeffi- 
cients f], (, and hence I will focus here exclusively on second order transport coefficients 
and review the existing results for relativistic field theories. Of particular interest are 
the values of the "relaxation times" r n , rn in the shear and bulk channel, respectively, 
since these have implications for numerical solutions of relativistic viscous fluid dynam- 
ics. 

For conformal fluids in 3+1 spacetime dimensions, r n is known for quantum field 
theories at weak coupling A < 1 [0, and for a particular quantum field example 
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(A/" = 4 SYM theory) at very large coupling A ^> 1 [p3fl , and for infinite coupling 



A ->■ oo 14, B51: 



lim t~ 

A^O 



5.9- 



7] 



lim r w 

A— >oo 



4 _21n2 + — C(3)A- 3 / 2 



e + P ' aToo "" V ' "' ' 8 ' J e + P 
For non-conformal fluids, Tn is known only in one particular example of a strongly 
coupled field theory |26] where 

i - ) . r. = i I -'21n:2i- // 



Remarkably, for these cases, the maximal propagation speeds ( |3.6|) turn out to be 
less than unity, indicating that the second order fluid dynamic equations of motion obey 
causality Furthermore, it was shown recently that causality in second order fluid 
dynamics is guaranteed by the causality of the underlying quantum field theory for the 



conformal field theory example of M = 2 SYM at strong coupling [27]. This seems to 
suggest that second order fluid dynamics may be a particularly useful approximation 
of quantum field theories in the appropriate limit. 



The parameter Ai is known for quantum field theories at weak coupling p3| , and 
for M = 4 SYM theory at very large, and infinite coupling O, p 



lim Ai 

A^O 



5.2- 



rj- 



lim Ai 

e + P a->oo V 4 ' v ' J e + P 

The parameter A2 is known for quantum field theories at weak coupling [23|, and for 
M = 4 SYM theory at infinite coupling [25] : 



2 + ^C(3)A- 3 / 2 



rj 



lim A2 

A^O 



lim A2 

A— >oo 



In 2^. 

7rT 



In attempts to calculate A3, the parameter has been found to vanish at weak and strong 
coupling [ 



25 



The parameter n has been calculated in weakly coupled SU (N) gauge theory p9| 
and for M = 4 SYM theory for large, and for infinite coupling fl4] , p4| : 

5s s / 145 



lim k 

A^O 



lim k 

8tt 2 T ' a^oo 4tt 2 T 
Note that k 7^ 2Ai, except for M = 4 SYM theory at infinite coupling. 



1 



-C(3)A 



-3/2 



For the strongly coupled theory example of Ref. [pq] , all second order transport 
coefficients for a particular class of non-conformal fluids may be extracted: 



£(1 



3c 



6 



3 



r* = -^(1-3^), A 4 = 
l-3c 2 ), 6 = #(l-3c 2 ), 



6 = f (1 - 3c 2 ) , £ 6 



& (1 - 3c *) 



(5.1) 
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As a non-trivial consistency check, one can proceed to evaluate the relations between 
£ 5 ,£ 6 and the other transport coefficient derived in section [Q| . Using T = 
k oc s/T, one finds 



and 



6 



(34 - 1: 



--M 



1 _ 6c 2 + 3^/) + 2T(1 - 3^)A 3 (2^ - B' 3 



-(1 

3c? 



3cf 



(« + A 3 <^) 



which precisely matches the values from Ref . [ 2B| . It would be interesting to check these 
relations for other examples of non-conformal field theories (e.g., by compactification 



of conformal field theories suggested in [30 



Furthermore, one can use the above results for the conformal fluid transport co- 
efficients for infinitely strongly coupled M = 4 SYM theory and evaluate the three 
relations Q,(gD,Q: 



4irr] 



A, 



7T 



Tf 



4:711] 



A, 



4tti] 



Bi 



1 

4 



(5.2) 



These values precisely correspond to those for the entropy current derived from the 
horizon of a black hole in AdS§ (see Ref. P21). This serves as another consistency- 
check of the above approach. 

It would be interesting to calculate the entropy current from the black hole horizon 
corresponding to the geometry considered in Ref. |26| and check that its form fulfills 



the conditions on Bi, B 3 , £? 4 , B 5 , B 6 , A 2 , A 3 , A 5 found in section 4.2 



Also, it could be possible to constrain the coefficient Ai for M 
existing results in a gradient expansion beyond second order . 



4 SYM using 



6. Discussion and Conclusions 

In this work, I have derived the structure of the energy-momentum tensor and non- 
equilibrium entropy current for a relativistic viscous fluid in curved space at vanishing 
charge density up to second order in gradients. I found that in this case there are 
15 possible transport coefficients multiplying second order gradient terms (10 of which 
vanish for conformal fluids), in addition to the bulk and shear viscosity coefficients 
arising at first order. Requiring the divergence of the (non-equilibrium) entropy current 
to be positive definite, two of these 15 coefficients are found not to be independent, 
and their relation to the other coefficients is specified. 

Also, the most general form of the entropy current out of equilibrium allows 13 
possible coefficient functions (8 of which vanish for conformal fluids) multiplying terms 
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of second order in gradients. Requiring positivity fixes 8 (conformal fluids: three) of 
these coefficient functions. 

These results resolve any ambiguities in the structure of the viscous fluid dynamic 
energy-momentum tensor and thereby clear the path for an extraction of the (second 
order) transport coefficients from quantum field theories. For instance, using Eq. ( |3.2| ) 
and Eq. ( |3.3| ) one can derive Kubo-like formulas for correlators in the fluid dynamics 
regime for non-conformal fluids. In addition, knowledge of the correlators and spectral 
functions to second order in gradients may be vital to extract first order gauge theory 
transport coefficients from lattice Quantum-Chromodynamics . 



Also, knowing the structure of the viscous fluid dynamic equations may be required 
to extend conformal fluid dynamics simulations of heavy-ion collisions J33], |34], |35|, |36| 
to reliably include effects from bulk viscosity. It could furthermore help to clarify the 
role of non-linear viscous damping of the r-mode instability in rotating neutron stars 
§§§ 

Moreover, the agreement ( |5.2| ) of the non-equilibrium entropy current from fluid 
dynamics and that derived from the horizon area of a black hole support the hope of 
linking these two concepts in more detail, which might have wide ranging consequences. 
In this context it is curious to note that for conformal fluids, where one parameter (A\) 
could not be fixed by requiring positivity of the entropy current alone, a diffeomorphism 
ambiguity prevents fixing of the same parameter on the gravity side. It is unknown 
whether A\ (and hence the fluid entropy current) cannot be fixed in principle in fluid 
dynamics, or whether this is only a manifestation of our lack of ingenuity (cf. the 
discussion in Ref. for the gravitational standpoint). However, I believe A\ should 
have a definite value, since nature must know the amount of entropy production when 
taking a system out of equilibrium. 

Furthermore, it seems that genuine non-equilibrium contributions to the entropy 
current are necessary to ensure consistency of viscous fluid dynamics with the require- 
ment that entropy may never decrease. However, it is unknown to me whether the pres- 
ence of these non-equilibrium contributions to the entropy current is an experimentally 
established fact. If not, relations such as Eq. ( |4. 8|) could potentially be used to attempt 
such an experimental verification, since it was found that the divergence of the non- 
equilibrium entropy current differs from the "naive" expectation V^S 1 ^ = ^fO ilv <J iW . 
E.g., this could be done by measuring the entropy of a conformal fluid in equilibrium, 
then subjecting the fluid to a perturbation a^ v (ideally one that is tuned such that 
(Da^) ~ — \<J^ U (V ■ u) to eliminate the dependence on the unknown coefficient Ai), 
switching off the perturbation and letting the system relax back to equilibrium, then 
measuring again its entropy. Setups reminiscent of this proposal are being pursued in 
cold atom experiments ETJ, but it is unclear to me whether these could be adapted to 
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test for non-equilibrium entropy effects. 

Finally, it should be possible to extend my analysis to include the case of non- 
vanishing charge density, allowing for the presence of charge /heat diffusion by account- 
ing for gradients in the chemical potential, as well as including structures built out of 
the Levi-Civita symbol that break parity This would be the last step in determining 
the equations of motion for a general one-component relativistic viscous fluid. 
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A. Conformal Transformations 

Here I collect the behavior of terms up to second order in gradients under Weyl rescal- 
ings (|2.7| ). For the Christoffel symbols one finds 

rj* -> T% - (8%w + 5 x d u w - g, u d x w) , 

from which the behavior of the covariant derivative of the fluid velocity can be calculated 
to be 

V^u v ->■ e~ w iy^u u + Und v w - g^Dw) . 

This in turn implies 

"'fiv ' & "'{lis ? 

V • u — y e w (V • u — 3Dw) , 

which can be used to derive 

Furthermore, using the Christoffel symbols one can calculate the transformation of 
the Riemann tensor to be 

R x -> R x - S x V„d a w + 5 x V^d u w + g^V x d a w - g^V x d u w , 
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so that 



u x u u R <jMT>v + V< M <9 CT> w , 
^ i?^ + 2V li d v w + g^ u V x d x w , 
R -> e 2w (i? + 6V A <9 A u?) . (A.l) 

Finally, one has 

V^lns^ V£lns + 3V>, 
and the above result for the Christoffel symbols implies 

V^V^ In s -> V^V^ In s + 3V^d u> w + 2V^ In sV^w + ^a, u Dw , 

In s ^ e 2w (V^ In s + 3V^w + 3 (V • «) - In sV^w) . 

B. Useful Identities 

This appendix contains some useful identities that were, e.g., used in calculating the 
divergence of the non-equilibrium entropy current. 
Time derivatives of various scalars and vectors: 

n^Dn^ = -t^n^ (v • «) - n^Wp , 

£>V£ln a = -Vi (V-u)- Q^ + fl^ V?lns+ (c 2 -^ (V-u)V^lns 

-^cM InsVl Ins, 
DV^V^ In s = -V£ (V • u) - u u R^V*i In s 

f 9 - 3r 2 \ 

-V^VI In s f <v + (V • u) (B.l) 

2rvHl _„. / * . l + 3 C 2-3(rfln C 2/rflns)\ 
+c 2 V^ In sVl In s ( <v - (V • u) 2 K - ^ >\ 

+ In S ( V • u) - \v{a„ - Vjn^ , 

-D [-R + 2« Q ^i2 aj8 ] = 2 (u„iO - RT Q<v + ^ (V • - 2c 2 s u u R^ In 

Identities involving curvature tensors and vorticity: 

= 2c 2 lV In s + (VjruJ - D ff (V • «) 
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V^Vf In s - In sVf In s ( c 2 - 1^ 



2c 2 V^V^ In a + \viV^ a - ± V£ (V ■ u) - 20^0* 



+ (V • u) 



1 



4 



-ty^"" (V • it) , 
4 - 3r 2 

__• fi^fi'"' (V • u) + 2% v Q? v o% . 



Identities for scalar Sr- 



\ (|5 7 # + |S 3 ^ - g|5 4 ^ + ^5 5 «" + 2c 2 <V (B.2) 

+3Vf + 6V£ - 6V£ + (4 - 6^)Vi* - 12cX - f V£) = 
(c 2 In sV^ In s + V^; In s) + ( V ■ u) + In s - § ( V ■ u) In s . 



C. V At 5n 0n _ eq for Non-Conformal Fluids 



The explicit divergence of Eq. (|4.9| ) is 



v /i u non-eq 



+^V^ (V ■ u) (2(1 - 3c 2 s )A 3 - B 3 ) + A 5 DR 
+V^ Q V>s {B 2 (B> 2 - c 2 ) + 5 3 (B 3 + c 2 ) - B 5 ) 

+^V£ (V ■ u) V£ In s (-4c 2 (l - 3c 2 ) A 3 - 6A 6 + 3A 7 (^ 7 - c 2 ) - 25 3 (^ - c 2 ) + 35 6 ) 



^ + RT (~ + A 3 ) + UaUpBT^P I ~ K T ^ + ^ + Br - 2A 3 



A 3 + 27/TV 



4 A 





2T 


ln s 






V T 


2Ai 




T 



+ A x — 2A 2 - 2A 3 + Bx- 2B, 



+c 2 V? Vf ln a - Ax - 2A 3 -B 3 + c- 2 B A 



+ Ax + Bx-2A 3 )+ a 



(V • u) (2r){r« 



12 



12ft + 3Cr n 



-Ai(l + 3A[) + A 3 (4 + QA' 3 ) - 6A 4 - Si (4 + 3B[)) 
+V< fl lnsVT Ins 



2c 2 (c 2 c^tv + A4 + c , (c , _ ^ + ^ _ Ag 



2T 
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+A 7 (A' 7 - c 2 s ) + clB'.B, + c 2 (c 2 -B' 3 - c' s )B 3 + (B' 4 - c 2 s )B 4 )] 
+ (V • u) [| (V • u) + i? (-| + (1 - A' 5 )A 5 + i(3A' 3 - 1)A 3 ) 

^-l(l-SA' 3 )A 3 -2A?j 
CTn " 36 + ^(1- 34)^3 + \(1-S^)A^ 



+u a upR aP 
+ (V • uf 



3T 

-^±^ - 1(1 - 6c 2 s + 3A' 2 )A 2 + f (2 - 3c 2 s - 3A' 3 )A 3 + 2A 4 
-(c 2 s -B' 2 )B 2 -±(l-9c 2 s )B 3 -B^ 
+V£ V£ In sc 2 (^ + 4(cl - c' s - A' 3 )A 3 - 2A A + (1 - c~ 2 A' 7 )A 7 + cfB^ 



In sV^ In s 



c 2 (c 2 - c> s )(t u + £ 4 , 2 



T 



+ -c*(c* - 2< - 3c 2 s A' 3 )A 3 



(C.l) 



+2c^ - c'JA 4 + 1(1 + 6c 2 s - 3A' 6 )A 6 + 1(1 + Sc^c* - A' 7 )A 7 - (<* - B' 6 )B 6 



where it is recalled that cl 



dins 



' * dins ' * dins 



din Bi 
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